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Abstract: The maintenance of a soil’s infiltration rate (IR) and field-saturated hydraulic conductivity 
(Ks) is crucial for the long-term sustainable functioning of wastewater-irrigated lands. However, an 
effective procedure for reliably measuring in situ soil Ks remains elusive. To address this issue, this 
study investigated the dualhead infiltrometer (DHI), a novel instrument for automatically 
determining IR and Kg, and compared it with a traditional double-ring infiltrometer (DRI) under 
various field conditions. In the initial phase, we optimized the procedure and settings for the DHIs 
in a cropland that was spray-irrigated with secondary-treated wastewater for decades in central 
Pennsylvania. Results showed that our optimized procedure, which used a single, long pressure 
cycle, yielded more robust measurements of IR than the originally recommended sequence of two 
short pressure cycles. The values of Ks measured by the DHIs with optimized settings were similar 
to those measured by DRIs under many (but not all) field conditions, due to their differences in 
infiltration surface areas, operational procedures, length of infiltration time, and soil spatiotemporal 
variability. Viscosity-corrected Ky on the irrigated cropland was 123.8 + 94.0 mm-h", higher than 
that on the adjacent non-irrigated cropland (103.2 + 94.6 mmh“), but the difference was not 
statistically significant, owing to the high degree of soil spatiotemporal variability and our limited 
number of measurements. Nevertheless, the higher Ks values measured on irrigated cropland reflect 
observed changes in soil structure (e.g., soil pore characteristics) that resulted from decades of 
irrigation. Seasonal variations in Ks values existed between winter and summer conditions, but IRs 
during all seasons remained much higher than the current spray-irrigation rate (4.25 mmh“), 
suggesting that the soil is still capable of handling the routine irrigation, even during winter. 
However, the coefficients of variation exceeded 67.0% across the field sites investigated, and the 
time periods covered by our measurements were limited. As this specific site is permitted to 
discharge treated wastewater year-round, caution must still be exercised to ensure that soil Kp 
remains high enough to prevent runoff generation, especially during winter frozen conditions. 


Keywords: hydropedology; soil physics; soil moisture; soil temperature; wastewater irrigation; 
infiltrometer 


1. Introduction 


Many soil physical, chemical, and biological factors influence infiltration rate (IR) and field- 
saturated soil hydraulic conductivity (Ks), both of which show significant variations across space and 
time [1,2]. Porosity and pore size distribution are among the major controlling factors of the rate of 
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water flow in soils [3]. For instance, greater macroporosity often gives rise to a higher IR [4]. High 
temperatures increase IRs by changing the viscosity and surface tension of the infiltrating water 
[3,5,6]. A higher antecedent soil moisture limits the IR [7,8]. Vegetation and land use can also alter 
soil characteristics, such as aggregation, organic matter content, compaction, and pore size 
distribution, all of which impact soil hydraulic conductivity [3,9]. For example, vegetation canopies 
and leaf litter can maintain the infiltration capacity of forest soils by protecting the soil surface from 
the impact of raindrops and preventing the formation of surface crusts, while roots induce the 
formation of macropores that improve the soil’s macroporosity and infiltration capacity [10,11]. 

In addition to spatial variations due to site-specific conditions, soil IR and Kp also experience 
temporal dynamics. For instance, hydraulic conductivity is often higher in the dry season than in the 
wet season [12,13], and higher in the summer than in the winter [14,15]. Soil infiltration at high 
latitudes or elevations is affected by frequent freeze-thaw cycles [16] and the presence of ice-rich 
layers, which can impede infiltration during winter [17]. 

Wastewater irrigation is a cost-effective and sustainable water resource management practice 
that is gaining attention as the world faces a looming freshwater crisis [18]. One priority of long-term 
wastewater irrigation operations is to maintain adequate soil infiltration capacity and Ks, which may 
vary with soil type, water quality, and management practices [19,20]. While some researchers 
identified benefits of applying treated wastewater irrigation on soil properties, such as reductions in 
bulk density, greater aggregate stability, and increases in organic matter content [20,21], other studies 
found no such improvements or identified adverse effects, such as reduced infiltration capacity [22— 
25]. In terms of soil hydraulic properties, both positive and negative influences were reported. For 
example, Hati et al. [20] and Vogeler [21] reported increases in saturated and unsaturated hydraulic 
conductivity, respectively, after treated wastewater irrigation, while Wang et al. [25], Coppola et al. 
[26], and Gharaibeh et al. [27] found that infiltration rates, permeability, and saturated hydraulic 
conductivity declined with the application of treated wastewater effluent, respectively. An increase 
in total porosity, a larger mean weight diameter of aggregates, and a lower bulk density could 
increase hydraulic conductivity, whereas increased retention parameters and bulk density, or pore 
occlusion due to migrated clayey material could reduce hydraulic conductivity [20,21,27,28]. The 
balance of these effects will determine the net change in soil hydraulic conductivity, which is likely 
to be highly dependent on local conditions and management practices. 

Since 1962, Penn State irrigated its secondary-treated wastewater on a specific area of cropped 
and forested land [29]. The naturally thick soils chosen for this purpose effectively acted as a “living 
filter” to provide tertiary—and final—treatment of the treated wastewater (Figure 1). During the 
period from 1962 through 1976, this land-based spray irrigation system was investigated as a small- 
scale test site. In 1983, it was expanded into a much larger and full-scale wastewater-treatment 
operation. As a result, Pennsylvania State University stopped discharging its effluent directly into 
the local watershed tributaries. Instead, the treated wastewater was spray-irrigated year-round to 
provide (1) tertiary effluent treatment before the treated wastewater enters into water bodies (water 
quality protection), (2) re-use of nutrients for crop production and tree growth (nutrient 
management), and (3) enhanced recharge to the groundwater to re-supply water to wells (water 
quantity conservation). 

Traditional wastewater disposal generally focuses on direct discharge and on treatment 
methods designed to achieve levels of effluent quality that are acceptable for discharge. However, 
regulations aimed at reducing water resource degradation require that environmentally sound and 
cost-effective non-discharge alternatives must be considered and used. Environmental consciousness 
and emerging technologies for non-discharge wastewater management alternatives evolved over 
time. Recycling wastewater for re-use and returning treated wastewater to groundwater aquifers by 
land application emerged as strategies to respond to this need. Non-discharge practices constitute 
good environmental management, and states in the United States (US) are required to adopt an anti- 
degradation policy that meets minimum requirements. They must include this policy as a required 
element of their surface water quality standards programs in order to gain federal approval of the 
standards. 
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In general, subsurface wastewater treatment applications offer numerous advantages [30]. These 
systems typically require a low degree of maintenance and the energy requirements are also low. 
Temperature equilibration of water is achieved during subsurface storage and excursions in water 
quality are buffered due to dispersion in the subsurface and dilution with native groundwater. The 
“living filter” allows reclaimed water to infiltrate slowly through the vadose zone, where sorption, 
filtration, and biodegradation can enhance the water quality (also called soil aquifer treatment [30]). 
To maintain soil infiltrability during operation with reclaimed water, and for mosquito control, 
wastewater spray fields are usually operated in alternate wet and dry cycles. As spray fields dry out, 
dissolved oxygen penetrates into the subsurface, thereby facilitating biochemical transformation 
processes, and organic material accumulated on the soil surface will decompose, allowing for the 
recovery of infiltration rates [31]. 

Following years of research in Pennsylvania’s northern climate, and experimentation with 
various application rates (2.54, 5.08, 10.16, and 15.24 cm/week), a 5.08-cm/week (i.e., two inches/week) 
rate was selected for year-round use in the “living filter” [32]. As a 5.08-cm/week spray rate is 
generally applied over a 12-h period, once per week (the rest of week, the soil is rested to dry out), 
which is equivalent to 4.25 mm of effluent per hour [33]. Each lateral at the “living filter” is only 
turned on for between 25 and 35 irrigations each year, resulting in an average annual effluent of about 
1500 mm [33]. Two large spray fields were carefully chosen for this purpose. These sites are underlain 
by thick residual soils derived from weathered carbonate rock and a deep water table [32]. Irregular 
surface topography, resulting from the differential weathering of carbonate bedrock, ensures 
retention, infiltration, and recharge of effluent. 

Many methods and infiltrometer designs were used to measure IR and to calculate Kg [34,35]. 
To alleviate some of the difficulties encountered with reliable and automatic measurements, 
dualhead infiltrometers (DHIs) [36] were developed based on falling-head measurements [37]. These 
DHIs measure Ks using two pressure heads, which are controlled by air pressure rather than the 
water ponding depth. To date, this automated device was subjected to limited field testing (e.g., 
Reference [38]). Previous studies typically applied the instrument using the default settings 
recommended by the manufacturer (e.g., References [39,40]). Therefore, its performance, optimal 
measurement procedures, and related parameterization under various field conditions need more 
field tests. 

The objectives of this study were to (1) compare the performance of the new DHI using original 
and optimized settings; (2) compare the results from the DHI with those from classic double-ring 
infiltrometers (DRI); (3) investigate the seasonal dynamics of Kg, especially under winter and summer 
conditions; and (4) determine if, and to what extent, the long-term irrigation of wastewater influenced 
Kg and its spatiotemporal dynamics, and particularly the long-term sustainability of the living filter. 


2. Materials and Methods 


2.1. Study Sites 


The study area is located in the Pennsylvania State University’s “living filter” in central 
Pennsylvania (Figure 1), at an elevation of about 335 m. The climate is a composite of the relatively 
dry midwestern continental climate and the more humid conditions found near the eastern seaboard 
[41]. Mean monthly air temperature peaks in July (22.2 °C) and bottoms in January (—2.8 °C), with 
daily minimum air temperature generally remaining below freezing from November through March. 
Mean annual precipitation is 932.4 mm, with the most precipitation occurring in May (101.6 mm) and 
the least in January (59.4 mm). Annually, the study sites receive a total of about 2.5 m of water 
(irrigation plus natural precipitation), resulting in a “northern rainforest” [33]. The living filter is 
situated on a thick (up to 100 m) residual limestone-derived Hagerstown silt loam soil [42]. It is well 
drained with moderate permeability [29]. The regional water table is roughly 60 m below the land 
surface [43]. 
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Figure 1. Penn State’s “living filter,” located in central Pennsylvania (a), with two fields—one is the 
State Game Land site (425 acres), and the other is the Astronomy site (182 acres); altogether, there are 
~66% forested and ~34% cropped lands. The zoomed-in map (b) shows the cropped portion of the 
Astronomy site and indicates the infiltration measurement locations (K1, K2, and K3). Also included 
are photos demonstrating the operation of the living filter during the (c) summer and (d) winter. 


Selected chemical properties of the treated wastewater are shown in Table 1 [44]. The no-till 
cropland is planted with an annual rotation of corn (Zea mays L.) and wheat (Triticum aestivum L.). 
Three sites were used to measure Kg (Figure 1): a control site, K1, located on the northwest margin of 
the cropland but outside the range of the irrigation; and two experimental sites, K2 and K3, located 
within the irrigated cropland. The distance between neighboring sites was approximately 100 m. 


Table 1. Average chemical characteristics of the municipal wastewater used for irrigation from 
November 2005 until August 2006 (n = 21) [44]. 


Parameter (Unit) Values 
Total alkalinity (mg-L“! as CaCOs) 179 
Biological oxygen demand (mg-L“) 8 
Specific conductance (umhos-cm*) 1109 
Total hardness (mg-L7! as CaCO») 245 


pH 7 
Total suspended solids (mg-L“) 5 
Nitrate (mg-L~! as N) 12 


2.2. Field-Saturated Soil Hydraulic Conductivity Measurements by DHIs 


2.2.1. Field Measurements 


The Kg was measured intermittently from February 2016 to June 2017 using DHIs (Figure 2a). 
The exact locations of infiltrations varied weekly to avoid repeated disturbance to the soil, but all 
infiltrations at each site were conducted within a 10-m2 area on the flattest topography available. 
Leaves, large vegetative litter, and snow were carefully removed during the preparation of the soil 
surface for infiltration measurement. The 5-cm-deep, 7.5-cm-radius insertion ring (Figure 2a) was 
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gently hammered into the soil to ensure good contact with the soil and minimal disturbance. It was 
then checked to ensure that it was level in all orthogonal directions. Prior to data collection, ground 
surface conditions were recorded, and antecedent soil volumetric water content (VWC) and soil 
temperature (Ts) were measured nearby at a depth of 10 cm using ECH20 STE sensors (Meter Group, 
WA, USA). To measure the temperature of the infiltrating water (Tw) during the experiment, a 
thermometer was secured inside the chamber of the DHI before the infiltrometer head was clamped 
onto the insertion ring. Air temperature (Tz) was taken from the Rock Springs site of the Soil Climate 
Analysis Network (SCAN), 16 km away from the living filter [45]. 

The default measurement parameters and the device program were tested and then compared 
with optimized settings during the first few winter months of 2016. Subsequent infiltration 
experiments adopted only the optimized parameters (Table 2). The default procedure (for wet, loamy 
sand) included a 15-min soaking period and two cycles of 15-min holding times at the high and low 
pressure heads [36]. Our optimized procedure was implemented to better ensure that a steady state 
of infiltration was reached (a prerequisite for reliable Ks measurements). It replaced the two short 
cycles with a single 35-min cycle at each pressure head (Table 2). The DHI recorded the water depth 
and pressure head (Figure 3), as well as the water flux per minute. During each infiltration 
experiment, the infiltrometer head was checked periodically to ensure that the seal was intact. If any 
sign of leakage or near-surface lateral flow occurred (for example, due to topography or frozen soil), 
the test was aborted. The ponding depth of the DHIs remained the same, at 5 cm, for both the lower 
and higher pressure heads, which was equal to the lower pressure head (5 cm). To achieve the higher 
pressure head (10 and 15 cm for the default and optimized setting, respectively), an air pump in the 
control unit (Figure 2a) was used to pump air into the sealed infiltrometer head (Figure 2a) to add air 
pressure. 


Figure 2. (a) Dualhead infiltrometers (DHIs) during a measurement (in triplicate) on the irrigated 


cropland, and (b) double-ring infiltrometer (DRI) during a measurement alongside a dualhead 
infiltrometer on the non-irrigated cropland. 
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Figure 3. Examples of the raw data recorded during measurements using a DHI with (a) the default 
parameters and (b) the optimized parameters, as given in Table 2. The “5 points adjacent-averaging” 
line smoothed the data by averaging consecutive data points within a moving window of five points. 


Table 2. Characteristics of eight infiltration processes with different measurement settings (default 


vs. optimized). CV_PH and CV_PL are the coefficients of variation (CV) for high and low pressure 


heads, respectively. CV_IH and CV_IL are the CVs for infiltration rates at high and low pressure 
heads, respectively. ISIp_H| and |Slp_L! are the absolute values of the linear regression slopes of the 
infiltration rate over time at high and low pressure heads, respectively. Different letters indicate a 


significant difference at a = 0.05. 


Parameters and Statistics Default (Original) Optimized 
Setting (Modified) Setting 
Soak time (min) 15 15 
Low pressure head (cm) 5 5 
arena High pressure head (cm) 10 15 
Hold time at each pressure head (min) 15 35 
Pressure cycles 2 1 
Total run time (min) 75 85 
CV_PH (%) 4,043.28 0.9+1.1? 
CV_PL (%) 1.1 +0.43 1.2+0.8° 
ES CV_IH (%) 7144.54 3.0+1.5° 
CV_IL (%) 6.043.94 4543.88 
|Slp_H! (mmh?) 75.5 + 42.38 50.9 + 61.98 
|Slp_L| (mm-h*) 29.2 + 33.7 2 17.6 + 18.4 è 
Number of test datasets used (see Figure 4) 8 8 


2.2.2. Quality Controls 


All data collected by the DHI were screened carefully to ensure quality. The coefficients of 
variation (CVs) of the pressure head and water level were calculated. As is common in statistical 
analysis, if a CV was greater than 15%, the data were considered unreliable and were subject to 
elimination [46]. In our case, this threshold was used to identify untrustworthy tests and 
measurements during which the pressure head, water level, or IR was not stable, possibly due to a 
poor seal. Considering that the IRs were recorded once per min, and this high frequency could record 
some high-frequency fluctuation of the IRs, some IRs may fluctuate within a range and did not show 
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a clear increasing or decreasing trend and, if they were resampled in 2 to 4 min, the CV of IRs could 
decrease greatly. Thus, if a measurement had a CV of raw IRs larger than 15% but lower than 50%, 
and the CV of resampled IRs was lower than 15%, then it was also thought to be reliable. Only the 
reliable measurements were used in further analysis. 


2.2.3. Calculation of Kis 


The Kg of soil was computed as follows [37]: 
Kiz (1) 


where i is the steady-state IR (mm-h-’), and F is a function that corrects for soil sorptivity and the 
geometric effects of the infiltrometer. The function F was calculated using the method of Reynolds 
and Elrick [47] as follows: 


A+D -2+2 


F=1+ $ 
Cd + C,b A 2) 


where D is the ponding depth (mm), d is the insertion depth of the infiltrometer (cm), b is the 
infiltrometer radius (mm), A is a shape factor which is equal to Cid + C2b (mm), C1 is 0.993, C2 is 0.578, 
and A is the macroscopic capillary length of the soil (mm). Since Kg was identical for both ponding 
depths (pressure heads), it was defined as 


iA iA 
Ks = = , (3) 
A+A+D, A+A+D, 
and the value of Kg can be solved as 
A(i, -i,) 
K, ==, 
f D -D, g 


where A is a constant for a given infiltrometer geometry. For the DHI, d = 50 mm and b = 75 mm; thus, 
A = 93 mm. As recommended by Decagon Devices [36], only the data from the last cycle were used 
in the Kg calculation, and the first two minutes of raw data after the pressure head changed were not 
used for calculating Ks. To ensure comparability, a similar subset of pressure and infiltration records 
at each pressure head was used to calculate Ks for the optimized measurements, even though the 
cycles took longer. 


2.2.4. Viscosity Correction of Ks 


All Kg values were viscosity-corrected to a standard temperature of 25 °C before comparing 
among different sites and different times to avoid the influence of changes in water effluent viscosity 
at different temperatures. Ks may be split into two factors [48] as follows: 


Kobe) 6) 


where k is the intrinsic permeability of soil (m?), and fis the fluidity of water (m~-s-). The f parameter 
is inversely proportional to viscosity, as 


a ae (6) 


where pw is the density of the water (kg:m™), g is gravitational acceleration (m-s~), 7 is the dynamic 
water viscosity (Pa's), and 1’ is the water’s kinematic viscosity (m%s). The corrected Ks was 
calculated as 


= K ie Nog (7) 
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where Ky_oz is an observed Kg value, and Kg_s7r is the Ks value corrected to the standard temperature. 


Mog and ex are the kinematic viscosity at the temperature when Kg_os was measured and at the 
standard temperature (25 °C), respectively, and each can be calculated as follows [49]: 


1825.85 


nr = 1.98404 x 1076 x e27347, (8) 


where 1. is the kinematic viscosity (m?-s~!) at temperature T (°C). 


2.3. Measurement of Ks by DRIs 


The DRIs used for this study had an inner ring radius of 3.02 cm and an outer ring radius of 5.40 
cm (Figure 2b). Each infiltration measurement was conducted for 85 min to match the total duration 
of the DHI procedure. The device was carefully pushed into the surface soil to minimize disturbance. 
After a soaking period, the water was manually poured into the device and the change in the water 
level was recorded every minute. The value of Kg was calculated according to Philip’s infiltration 
equation [50], 


i = St + At, (9) 


where i(t) is the cumulative infiltration at time t, S is the sorptivity of the soil, and A is a hydraulic 
conductivity parameter. Previous studies [51,52] suggested that the value of A is between 2/3 Kg and 
1/3 Kg, but usually closer to 1/3 Kg. For this study, the Ks measured by the DRI (Kp_pr) was calculated 
as Kys pr = 3A. 


2.4. Statistical Analysis 


The mean IRs and Kg values of a period or a site were calculated along with their standard 
deviations and presented as means + standard deviation. Independent sample t-tests with a set to 
0.05 were used to compare the Tu, Ts, VWC, IR, and Ks values measured at different sites and times, 
as well as the CVs of pressure heads and infiltration rates, and the absolute linear regression slopes 
of the infiltration rate over time. To satisfy the assumption of normality, log-transformed 
corresponding data were used in the t-tests. This analysis was completed using SPSS 17.0 (SPSS Inc., 
Chicago, IL, USA). Linear regressions between Ks measured by DHIs and DRIs, between Kg and Tu, 
and between VWC and Kgsrp were fitted using OriginPro 2019b (OriginLab Corporation, 
Northampton, MA, USA.). 


3. Results and Discussion 


3.1. The Influence of DHI Holding Time 


The default (original) and modified (optimized) DHI measurement parameters resulted in 
different infiltration patterns (Figure 3). Pressure head was not steady during the short high-head 
cycles, even during the second cycle (Figure 3a), and the CV of the high-pressure head was 
significantly higher for the default procedure than for the optimized procedure (Figure 3b, Table 2). 
Since higher pressures led to higher IRs, the performance improvement of the optimized procedure 
was more pronounced during the high-head cycles (Table 2). In addition, since the high pressure 
head was achieved by increasing air pressure rather than by pooling additional water [36], the quality 
of these tests depended on complete airtightness of the infiltrometer head. By contrast, the pressures 
and IRs during the low-pressure periods were not significantly different between the default and 
optimized programs (Table 2). This pressure was set equal to the nominal ponding depth; thus, it 
was achieved with actual water-level head, rather than extra air pressure. The duration of the default 
procedure afforded 15 measurements of IR, but the first two measurements at each pressure head 
were not used in the Kg calculations, as recommended by the device manual. However, the third and 
fourth records were also often unsteady as they continued to be affected by the transition between 
heads (Figure 4). Moreover, IR could be expected to fluctuate in response to changes in the physical 
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conditions of the soil, such as the melting of ice, the release of trapped air bubbles, or the 
breakthrough of water into connected macropores. These processes could distort the data for several 
minutes at a time, making the true steady-state IR difficult to determine with such a short available 
time window. Anomalous time-steps were excluded from Kg calculations, which sometimes left 
fewer than 10 records from which to calculate Kg. 

Figure 5 contains time-series plots of all successful infiltrations collected with the optimized 
procedure. These results are grouped by season and antecedent VWC conditions. Winter data were 
collected from November through March, and summer data were gathered in June and July. The soil 
was considered wet when its antecedent VWC was higher than 0.25 m3-m-3, while values below this 
cut-off were considered dry. Several anomalous data points are visible in these plots (e.g., 
measurement #22 at 39 min and measurement #33 at 73 min; Figure 5c,d). 
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Figure 4. Infiltration rates of experiments conducted in winter months with (a) the default settings 

and (b) the optimized measurement parameters, as given in Table 2. 


The optimized procedure, with its single long cycle, yielded better infiltration results than the 
default procedure for the silt loam in the living filter (Figure 2b). A long cycle increased the likelihood 
that the test reached a steady pressure head and IR, and it provided more of a buffer against 
fluctuations caused by transient changes in soil conditions (Figure 2b). Field applications of DHIs in 
soils heavier than a sandy loam may benefit from this improved approach: a 15-min soak time 
followed by 35 min held at a high pressure head, and then 35 min held at a low pressure head (Figure 
3b). The IRs measured during the final 15 to 20 min of the hold time at each pressure head were 
generally usable for the calculation of Kg. 
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Figure 5. Infiltration rates under (a) wet and (b) dry soil moisture conditions and in winter (November 
to March), and under (c) wet and (d) dry soil moisture conditions in summer (June and July), and 
under (e) special conditions in winter (with snow cover, Ts < 0 °C, or Ta < 0 °C). Wet soil refers to an 
antecedent soil volumetric water content (VWC) higher than 0.25 m3-m°3, while dry soil had a VWC 
below this value. T4 is air temperature, and Ts is soil temperature. Each line represents a single 
infiltration run, with the legend numbered values corresponding to the “Number” column in Table 
3: 


3.2. Comparison of Ks Measured by DHIs and DRIs 


The values of Kg measured by DHIs (Kg_px) and DRIs (Kg_pr) on the same day are compared in 
Figure 6. Nine of 14 Kg results were similar between the two methods (distributed along the 1:1 line). 
However, the other five Ks_p1 values were much higher than Kg_pr. The regression line of all 
measurements has a slope of 1.15 mm-h" and an intercept of 13.45 mm-h"!, showing that, on the 
whole, the DHI measured higher Ks values than the DRI. 

This difference between the performances of the two instruments may be explained by the larger 
footprint of the DHI. Kumar [3] also found that measured IR values were lower when using smaller 
DRIs. Lai and Ren [53] observed that smaller inner rings resulted in greater variability in Kg 
measurements. The infiltration area of the DHI is 13.8 times that of the inner ring area of the DRI used 
in this study. This makes it less sensitive to the spatial heterogeneity of the soil and probably closer 
to the representative elementary volume of the experimental site, which reduces the effects of the 
measurement scale [53,54]. 
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Figure 6. Relationship between field-saturated hydraulic conductivities measured concurrently by 
the DHI (Ks_pu) and DRI (Kg_pr). The solid dots indicate divergent measurements between the two 
methods, while the open dots, distributed along the 1:1 line (dashed), indicate similar results. All data 
were included in the regression line (solid). The equation, coefficient of determination (R?), and the 
significance value (p) are also shown for the all-data regression (solid line). 


On the other hand, smaller rings require a smaller area of flat topography and are more easily 
operated in sloped areas. The DHI uses a peristaltic pump in the control unit and a water-level sensor 
in the infiltrometer head (Figure 2a) to maintain the water level and measure the water flux. This is 
operator-friendly, but it also relies on an air pump to control the pressure head, which adds to the 
complexity and uncertainty of taking measurements of Ks. During the collection of this dataset, about 
40% of all measurement attempts failed, usually as a result of a compromised seal between the ring 
and the infiltrometer head, or a CV of the IR higher than 30%, leading to the inconsistent numbers of 


repetitions listed in Table 3. 
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Table 3. Date, site, air temperature (T4), soil temperature (Ts), antecedent soil volumetric water content 
(VWC), infiltration rates at the higher and lower pressure head (IR_H and IR_L, respectively), and 
field-saturated soil hydraulic conductivity (Ks) of measurements using the optimized parameters 
defined in Table 2. Corresponding raw infiltration data are graphically presented in Figure 5, and 
different backgrounds from top to bottom in this table correspond to Figure 5a, b, c, and d, 
respectively. 


Effective Ta Ts VWC IRH IR_LL Kfs 


Number Measurement Date Site # 


Repetition °C °C ë mèm? mmh! mmh! mmh 


32 3 June 2016 K2 1 Vl 22,9) 0.18 311.4 680.1 304.9 
33 3 June 2016 K2 2 6.0 22.2 0.18 284.0 643.0 325.0 
34 12 June 2017 K2 1 26.7 26.2 0.18 13.2 45.0 40.6 

35 12 June 2017 K2 2 26.7 26.2 0.18 45.4 86.6 48.7 

36 12 June 2017 K3 1 26.7 26.0 0.18 72.0 184.3 116.1 
37 3 June 2016 K1 1 Zee 22.9) 0.23 284.0 643.0 325.0 
38 3 June 2016 K1 2 4.6 20.5 0.23 367.2 730.5 284.4 


3.3. Seasonal Dynamics of IRs and Kp, and Their Influencing Factors 


Under most conditions, IRs decreased over time during the 15-min soaking period (Figures 
5a,c,d), with the notable exception of dry, winter conditions, when most measurements indicated a 
gradually increasing IR (Figure 5b). The increasing IR during the soaking step in cold, dry soil could 
be explained by an initial period of hydrophobicity caused by low temperatures and dry soil 
aggregates [55,56]. This effect would become less prominent as soil moisture rises during the soaking 
period. Under high pressure heads, the behavior of the IR depended on antecedent soil wetness: 
when the soil was wet, infiltration was likely to decline over time; when the soil was dry, IR tended 
to increase. 

During winter, when the soil was wet, the IR values clustered into two groups (Figure 5a). One 
group had higher IRs of 236.5 + 34.6 and 126.5 + 23.7 mmh“ at the high and the low pressure heads, 
respectively. The other group had lower IRs of 83.0 + 21.1 and 49.2 + 15.4 mm-h* at the high and the 
low pressure heads, respectively. These pronounced differences may have been caused by spatial 
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heterogeneity of soil properties (e.g., macroporosity), as the exact location of the experiment changed 
weekly. Also, it could be attributable to unidentified systematic errors. As reported by Gamie and De 
Smedt [57] and many others, great variation in hydraulic conductivity, spanning four orders of 
magnitude, was observed within an area of 120 m x 120 m. They attributed this to the combination of 
sampling and measurement errors, randomness, preferential flow, and soil heterogeneity. 

When the soil was dry and cold, we measured IRs that varied from 37.3 to 298.3 mm-h" at the 
high pressure head and from 16.1 to 142.3 mm-h" at the low pressure head (Figure 5b). During 
summer, when the soil was wet, the IRs at high and low pressure heads were 160.3 + 119.4 and 74.7 
+ 61.3 mm-h", respectively (Figure 5c). When the soil was dry during summer, the IRs at high and 
low pressure heads were 423.0 + 286.2 and 207.1 + 122.1 mm-h“, respectively (Figure 5d). Three of the 
measurements under dry summer conditions were much higher than the other measurements. 

During special winter conditions, such as periods with snow cover, air temperature (Ta) slightly 
lower than 0 °C (-2.6 °C < T, < 0 °C), and soil temperature (Ts) slightly lower than 0 °C (-0.6 °C < Ts < 
0 °C), IR values were 181.1 + 95.5, 193.3 + 87.6, and 207.7 + 79.8 mmh at the high pressure head and 
99.2 + 54.4, 96.7 + 44.3, and 109.4 + 37.9 mm-h"” at the low pressure head, respectively (Figure 5e). 
These [Rs were not significantly different from results obtained on other experimental days during 
winter (Table 4). A measurement on 3 February 2017 detected the highest IR during winter, despite 
the presence of snow cover, as well as cold and dry conditions (T: = -0.05 °C, Ts = -0.3 °C, and VWC 
= 0.20 mêm; Figure 5b,c). These results indicate that low air temperature, slightly frozen soil, and 
the presence of snow cover do not seem to impede infiltration at the study site. 
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Table 4. Characteristics of infiltration processes under different conditions. Ta is air temperature, Ts is soil temperature, and VWC is initial soil volumetric water 
content. IR_H and IR_L are the infiltration rates measured at high and low pressure heads, respectively; Ks_os and Kp_str are the initial observed and viscosity- 
corrected field-saturated hydraulic conductivities, respectively. Different letters in the same column indicate significance groupings at a = 0.05. 


Toe . Ta Ts VWC IR_H IRL Kfs_oB Kfs_stP 
Conditions Number of Times °C °C E RE anha mmha eae 
Winter, dry soil 6 4.74494 154194 0.2140.022 133.6 + 94.3 ə 56.4 + 45.8 4 70.7 + 46.3 a 88.3 + 60.3 4 
Winter, wet soil 9 2.7 £3.34 2.141.372 0.3140.01> 151.2+852 83.6 + 44.7 @ 62.2 + 45.0 88.0 + 57.6 
Summer, dry soil 7 17.1410.5% 23.8+23> 0.1940.032 423.0+286.2> 207.14122.1> 184.84141.0> 206.8+147.7> 
Summer, wet soil 17 18.942.7> 19.241.5° 0.30+0.01> 160.34119.42 74.7+61.34 77.5 + 57.4 a 90.1 + 63.9 a 
With snow 6 1.74532 0.24044 0.27+0.06° 181.1495.5% 99.2 +54.4 ab 75.1 + 45.2 a 111.8 + 62.8 a 
Ta< 0°C 4 -1.7+1.2a 014047 0.2840.05> 193.3+87.6% 96.7 +44.3 əÞ 88.5 + 47.33 123.3 + 63.3 əÞ 
Ts< 0°C 4 -0.14342 -0.2+0.34 0.28+0.06° 207.7+79.8%æ 109.4+37.9% 88.5 +47.4a 135.6 +60.9 ə 
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During winter, the mean Kg_oz was 70.7 + 46.3 mmh when the soil was dry and 62.2 + 45.0 
mm-h when it was wet. The corresponding summer values were 184.8 + 141.0 mm-h" and 77.5 + 
57.4mm-h", respectively (Table 4). After the viscosity correction, the mean winter Kp_ste values were 
88.3 + 60.3 and 88.0 + 57.0 mmh when the soil was wet and dry, respectively, while the 
corresponding summer values were 206.8 + 147.7 and 90.1 + 63.9 mm-h=, respectively (Table 4). In 
winter with snow cover, Ta < 0 °C, and T; <0 °C, the mean Kg_os values were 75.1 + 45.2, 88.5 + 47.3, 
and 88.5 + 47.4 mm-h“, respectively, and the corresponding values of K_srp were 111.8 + 62.8, 123.3 + 
63.3, and 135.6 + 60.9 mm-h, respectively (Table 4). Independent sample t-tests showed that the Kg_oz 
and Kg_sre values during dry summer conditions were significantly higher than the others (Table 4). 

The living filter is permitted by the Pennsylvania Department of Environmental Protection 
(PADEP) to operate year-round, rather than using storage during the non-growing season [33]. 
However, frozen soils are an important management concern in the living filter, as irrigation onto 
frozen soils could lead to ponding and surface runoff; the latter conditions could potentially lead to 
nutrient-rich water entering nearby surface waterbodies. To address this management concern, we 
collected real-time monitoring data since 2008, which indicated that frozen soil conditions do not 
seem to noticeably impede deep infiltration or significantly alter vertical flow. This implies that 
irrigation during frozen soil conditions may not adversely affect effluent treatment [58,59]. This is 
probably attributable to warmer temperatures associated with the effluent when the effluent was 
sprayed on frozen soils intensively and continuously, leading to locally permeable frozen surfaces. 
The depth and temporal extent of frozen conditions, as well as the infiltration status and subsurface 
water fluxes under frozen conditions, however, require further investigation. 

There were significant linear relationships between the observed Kp and Tw at both the irrigated 
and non-irrigated sites (e.g., K2 and K1, respectively; Figure 7). However, the viscosity-corrected Ky 
showed no significant relationship with Tw, suggesting that the temperature correction was effective. 
Additionally, no significant relationships were found between Kg_oz, Ks_sre, and soil temperature at 
K1 or K2. 


600a) K2 = —— K, op748.31+5.53T „ R= 0.10, p=0.03 

5004 o m Kp srp=114.62+2.87T „ R= 0.00, p=0.32 
~~ 400 j SiR 
a o 
: 300 | g i ê a 
~g 200 ] H nee a eg sae cB ge —a 
~ 100] OS a è 

0 $ g= 8 a 8 a a k a z 
5 10 15 20 25 

600 (b) KI = —— K, oy=-7.04+10.72T,, R°= 0.34, p=0.008 
ao g asss Kg srp=46.50+8.817,,, R= 0.34, p=0.061 
~~ 400 Aa o 
a o 

300 | 3 8 4 
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1004 — z fe 8 
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Figure 7. The relationship between Kg and infiltrating water temperature at (a) K2 and (b) K1 sites. 
Kg: field-saturated hydraulic conductivity; Ks_os: observed Kg; Kg_srp: Kg that was viscosity-corrected 
to the standard temperature of 25 °C; Tw: water temperature. 


The observed relationship between Ks and Tw is consistent with previous studies [6,15,60,61], but 
the relationship between Kg and soil temperature is different from the findings of Clancy and Alba 
[48], who found a significant linear relationship between the two variables. They suggested that soil 
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temperature could affect air viscosity and permeability, as the pore spaces of the field saturated soil 
are still between 5% and 20% occupied by air [62,63]. Air permeability is related to air viscosity and, 
thus, hydraulic conductivity correlates with soil temperature. However, we did not observe this 
association at the study site. Levy et al. [61] demonstrated that the effects of temperature on hydraulic 
conductivity were influenced by soil types. The results of this study suggest that the differences in 
Kgs at different temperatures were mainly controlled by the kinematic viscosity of water. 

We observed a significantly negative linear relationship between Kg_srp and antecedent VWC 
during winter, and across the whole experimental period, while such linear relationships were 
insignificant during summer (Figure 8). Similar linear relationships between Ks and antecedent VWC 
were found in many studies. Because wet soils have reduced soil sorptivity, subsoil drainage is often 
restricted [13,64]. In addition, high soil VWC tends to cause clay-rich soils to swell, constricting soil 
cracks, decreasing soil macroporosity, and reducing macropore flow [4,13]. 
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<~ 400 j 400 ° 3 
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Figure 8. Relationships between antecedent soil volumetric water content (VWC) and viscosity- 
corrected field-saturated hydraulic conductivity (Ks_str) at K2 during the entire study duration (a), in 
the winter (b), and in the summer (c). The shaded region in (c) encompasses notably low Ks values 
with the low soil VWC in the summer (included in the regression). 


Two Kj measurements on 12 June 2017 were extremely low. The surface soil was very dry, 
although the recorded VWC at 10 cm was at 0.18 mèm, which is not especially low. A similar result 
was observed on 5 July 2016, when the lowest soil moisture in the experimental period was recorded, 
at 0.14 m3-m°, and a relatively low Kg value was observed (Figure 8c). This reduced water infiltration 
may be attributable to water repellency of the soil [65,66], but it is also possible that an antecedent 
heavy storm or spray irrigation that fell on freshly harvested land during summer formed soil crusts 
which reduced the infiltration capacity of the surface soil [67]. 

The soil at the living filter tends to be cold and wet during winter [68], and this combined effect 
could explain the lower Kg values in winter than in summer. Similar seasonal variations of Ks were 
also observed in other regions and should be considered when planning wastewater treatment 
management and best management practices [12,14,15]. 


3.4. Long-Term Effects of Wastewater Irrigation on Kp 


In the irrigated cropland, at site K2, the observed minimum Kg was 13.2 mm-h and the 
maximum Kg was 311.4 mm-h"!, with a mean + standard deviation (SD) of 97.9 + 78.4 mm-h"! and a 
CV of 80.1%. There were only six valid measurements at K3, with a minimum and maximum Kg of 
37.4 and 174.3 mm-h-1, respectively, and a CV of 61.0%. Measurements at K3 were carried out in June 
2017 and had a mean Kg of 95.0 + 57.9 mm-h-, without a significant difference from the results of K2 
during the same time period. In the non-irrigated cropland, K1, the minimum and maximum Kg 
values were 14.8 and 367.2 mm-h“, respectively, with a mean + SD of 85.1 + 87.4 mm-h" and a CV of 
102.7%. The minimum of all results at all three sites was higher than the nominal application rate of 
the treated wastewater irrigation (4.25 mm-h-’). 

All viscosity-corrected Kg_stp values from K1 and K2 were compared, and no significant 
difference (p = 0.357) was found (Figure 9). At site K1, the mean and SD of Kg_stp was 103.2 + 94.6, 
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while the results at site K2 were slightly higher at 123.8 + 94.0 mm-h". During winter, Kg_stp values at 
K1 and K2 were 80.4 + 52.9 and 107.0 + 76.0 mm-h“1, respectively, without a significant difference (p = 
0.442). During summer, Kg_stp at K1 and K2 were 127.1 + 133.9 and 122.5 + 101.7 mm-h", respectively, 
also without a significant difference (p = 0.737). 


= Winter dry ¢ Winterwet 4 Summerdry vy Summerwet ® Others 


400 p=0.357 
: maximum 
| Density plot 
= 300 \ /|\Width = frequency @ a 
4 1 150% 
S 200 
al 
xX 
100 
0 


K1 K2 


Figure 9. Viscosity-corrected field-saturated hydraulic conductivities (Ks_str) in the cropland without 
(K1) and with (K2) irrigation, including a total of 23 and 29 valid measurements, respectively, during 
the experimental period. 


These results appear to be a bit different from previous findings in the same area. Sopper and 
Richenderfer [41] measured Ks using DRI and found that there was a significant increase in 
infiltration on the irrigated corn-cropped land. Walker [69] also reported a higher hydraulic 
conductivity in irrigated areas than in control areas, as measured by the constant head method in a 
laboratory setting. Larson [70] measured hydraulic conductivity using a tension infiltrometer and 
showed that hydraulic conductivity was relatively higher in irrigated lands than in non-irrigated 
lands, especially at 3 and 6 cm of tension, suggesting that the increased hydraulic conductivity was 
mainly attributable to an increase in the abundance of soil pores with radii of 0.025 to 0.050 mm. Due 
to the difficulty of data acquisition, the high failure rate of the tests, and the inherent spatiotemporal 
variability of the measured parameter, this study was too statistically underpowered to definitively 
answer the question of the effects of long-term irrigation on Ks. However, especially in light of other 
data collected at the same locations, it is likely that the increases observed here, although statistically 
insignificant, reflect real conditions in the soil. The insignificant results were due to limited replicated 
measurements at limited sites, combined with the considerable spatiotemporal variability of the 
measured parameter. 

Several explanations were proposed to explain increased hydraulic conductivities in irrigated 
soils. Firstly, the nutrient-rich treated wastewater encourages plant growth and soil organic matter 
accumulation. A lush crop canopy offers protection from raindrop impacts, thereby decreasing 
surface soil compaction or crusting, while the activity of plant roots improves soil structure and 
increases macropores. Macropores may contribute nearly 85% of the variation of hydraulic 
conductivity or total infiltration [71,72], as suggested by the dramatic increase in the incidence of 
preferential flow at irrigated sites in the living filter [42]. In addition, the increase in nutrients from 
the effluent would tend to improve soil physicochemical properties. Treated wastewater irrigation 
was found to increase organic matter content and soil pH [29], while lower bulk densities and higher 
soil pH values and electrical conductivities were measured on irrigated croplands, as compared to 


Water 2019, 11, 1632 18 of 22 


non-irrigated croplands [21,70]. These positive effects could result from higher organic matter inputs, 
increased soil fauna abundance, and plant root growth, as well as greater freeze-thaw cycling [73- 
75]. Many lines of evidence suggest that decades of spray irrigation caused physical adaptations in 
the irrigated soils at the living filter, which can accommodate much large volumes of water (beyond 
natural precipitation). 


4. Summary and Conclusions 


In situ measurements of two critical soil hydrologic parameters (IR and Kg) at the living filter 
wastewater tertiary treatment system at the Pennsylvania State University were conducted during 
both winter and summer months using newly designed dualhead infiltrometers (DHIs) following 
both original (default) and optimized (modified) operation settings. The results were also compared 
with traditional double-ring infiltrometers (DRIs) measurements under various field conditions. The 
results of this study support the following conclusions: 


(1) For the silt loam soil in our study site, the optimized DHI procedure, with one long cycle, 
achieved steadier pressure heads and infiltration rates and, therefore, yielded more credible Kgs 
values than the original two short-cycle default setting. 

(2) The values of Ks measured by the DHI were, for the most part, similar to those measured by 
DRI, but were higher under some conditions. This is likely due to the smaller infiltration area of 
the DRI we used, which possibly missed soil macropores. 

(3) The values of viscosity-corrected Kp from irrigated croplands were higher than those from non- 
irrigated croplands, but the difference was statistically insignificant (p = 0.357) due to high 
spatiotemporal variability and small sample sizes. Nevertheless, these results do suggest that 
the long-term application of wastewater irrigation at the site enhanced Ks. Measured IR and Kg 
values were all much higher than the actual spray irrigation rate under all conditions (4.25 
mm-h~). From this perspective, the current application of treated wastewater irrigation does not 
seem to be at high risk of causing overland flow, even under winter conditions with frozen 
ground. 

(4) Values of IR and Ks were lower during winter than during summer, as soil conditions tend to 
be cold and wet in the winter and warm and dry in the summer. Infiltration processes during 
winter should, therefore, receive extra attention when implementing wastewater irrigation 
systems, such as the one investigated here at Penn State’s living filter. 

(5) Values of IR and Kgs showed great spatiotemporal variability and some uncertainty, but they 
were clearly influenced by antecedent soil moisture and temperature conditions. Surprisingly, 
snow cover, cold air, and frozen soil did not show particularly significant impacts on winter IR 
and Ky in this study. This is probably attributable to warmer soil temperatures associated with 
the thick layers of effluent that tend to build up on frozen soils after long periods of intense 
irrigation. This aids in preventing runoff in the living filter system by encouraging locally 
permeable frozen surfaces. The depth and temporal extent of frozen conditions, as well as the 
infiltration status and subsurface water fluxes under frozen conditions, require further 
investigation. 
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